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We report complementary experimental, numerical and theoretical study of turbulent coflow, 
counterflow and pure superflow of superfluid '*He in a channel, resulting in a physically transparent 
and relatively simple model of decaying quantum turbulence that accounts for interactions of coex¬ 
isting quantum and classical components of turbulent superfluid ^He. We further offer an analytical 
theory of the energy spectra of steady-state quantum turbulence in the counterflow and pure super¬ 
flow, based on algebraic approximation for the energy fluxes over scales. The resulting spectra are 
not of the classic Kolmogorov form, but strongly suppressed by the mutual friction, leading to the 
energy dissipation at all scales, enhanced by the counterflow-induced decoupling of the normal- and 
superfluid velocity fluctuations. 

PACS numbers: 67.25.dg, 67.25.dk, 67.25.dm 


I. INTRODUCTION 

Flows of quantum fluids - such as ^He below T> ~ 
2.17 K - displaying superfluidity and the two-fluid be¬ 
havior, offer a challenging held of fundamental research 
combining quantum physics and fluid dynamics^^^. The 
phenomenological two-fluid model, suggested by Landau 
and Tisza, describes dynamics of superfluid ^He in terms 
of interpenetrating normal and superfluid components 
that have their own densities, Pn{T), Ps{T), and velocity 
fields, Un{r,t), Us{r,t). In this paper, we consider finite 
temperature - above about IK- where the normal com¬ 
ponent behaves as a classical fluid with the kinematic 
viscosity while the superfluid component is invis- 

cid, r's = 0. 

Due to the quantum mechanical restriction, the cir¬ 
culation around the superfluid vortices is quantized in 
integer values of k = h/m ~ 10“^ cm^/s, where h is the 
Plank constant and m denotes the mass of ^He atom. 
The singly quantized vortices usually arrange themselves 
in a tangle that can be characterized by vortex line den¬ 
sity (VLD) £, i.e., total length of the quantized vortex 
line in a unit volume. The dynamical behavior of the 
tangle constitutes an essential ingredient of quantum tur¬ 
bulence, the turbulence occurring in quantum fluids dis¬ 
playing superfluidity. The quantization of circulation in 
the superfluid component results in appearance of char¬ 
acteristic “quantum” length scale: the mean separation 
between vortex lines, £ = which is typically or¬ 

ders of magnitude smaller than the scale A of the largest 
(energy containing) eddies^’®. 

There is a growing consensus that the quantiza¬ 
tion of vortex lines can be neglected at large scales £ 
and that quantum turbulence is similar to classical turbu¬ 


lence if excited similarly, for example, in a steady chan¬ 
nel flow by a pressure drop®“^^ or when decaying be¬ 
hind a grid^^’^^. The reason is that the interaction of 
the normal-fluid component with the quantized vortex 
tangle leads to a mutual friction force^’^’^^ “which cou¬ 
ples together Un{r,t) and Us{r,t) so strongly that they 
move as one fluid” . On the other hand, at small length 
scales R "^ £, the quantization of vortex lines cannot be 
neglected and turbulence in superfluids has essentially 
quantum character. 

The pipe and channel flows of viscous fluids belong to 
the class of most extensively studied classical flows^®’^^. 
As for pipe and channel flows of superfluids, by combin¬ 
ing mechanical and thermal drive, a rich variety of two- 
fluid turbulent flows with different direction and flow ra¬ 
tio of the two components may be generated, representing 
a very complex superfluid hydrodynamics system^^d8.i9^ 
see Fig. 1. The classical-like mechanical forcing (e.g., by 
compressing a bellows) results in a coflow, the closest 
analogue to classical viscous channel flow, where both 
components move, on average, with the same velocity: 
the mean normal-fluid velocity is equal to the mean 
superfluid velocity Us, Fig. 1, right. However, due to 
quantum-mechanical constraints on the superfiow, gener¬ 
ated by individual lines, the velocity fields are completely 
different from the normal-fluid motions at scales Rf^£. 

The normal and superfluid components of ^He may 
also be made to flow relative to each other with a non¬ 
zero counterflow velocity 

Uns = Un-Us^0 . ( 1 ) 

The thermal counterflow, first systematically investigated 
in pioneering experiments by Vinen^*^, is easily generated 
in a channel with one of its ends sealed and equipped with 
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FIG. 1: (color online) Flow channels for the study of counter- 
flow, pure superflow, and cofiow. S and N stand for superfluid 
and normal components. Counterfiow is produced thermally 
by a heater. Superflow and coflow are driven mechanically 
by a bellows. The turbulence is probed in the middle of the 
channel by second-sound, excited and detected by mechanical 
vibration of a porous membrane. 

a heater and open at the other end to a superfluid he¬ 
lium bath (see Fig. 1, left). Here both components move 
relative to the channel walls. The heat flux is carried 
away from the heater by the normal fluid alone, and, by 
conservation of mass, a superfluid current arises in the 
opposite direction: 

Pn Un + PsUs = 0 ■ (2) 

In this way the counterflow velocity Uns, proportional to 
the applied heat flux, is created along the channel, soon 
accompanied by a tangle of vortex lines. In pure super¬ 
flow, sketched in Fig. I, middle, superleaks (i.e., filters 
located at the channel end with sub-micron-sized holes 
permeable only to the inviscid superfluid component) al¬ 
low a net flow of the superfluid component in the channel: 
C/g 0, while the normal-fluid component is remaining, 
on average, at rest: fTn = 0. In both cases the fields 
Us{r,t) and u^{r,t) are expected to be different at all 
scales. Thermal counterflow and pure superflow there¬ 
fore represent two special cases of counterflow, charac¬ 
terized by non-zero difference in mean flow velocities of 
the superfluid and normal components. 

In this paper, we report on complementary experimen¬ 
tal, analytical and numerical studies of the VLD decay, 
C{t), of three categories of turbulent channel flows of su¬ 
perfluid ^He, aiming to characterize the quantitative dif¬ 
ference between the statistics of turbulence in the coflow 
on one hand, and the counterflow and superflow on the 
other. Our study is based on detailed analysis of the VLD 
decay from three different initial values of Cq ~ 10^, 10® 
and 10® cm“^, in all three types of the flows and obtained 
at different temperatures^^ and represents therefore sub¬ 
stantial extension of recent Letter by Gao et , which 
considered the decaying counterflow. The experimental 
technique is shortly reviewed in Sec. IIA and sketched in 


Fig. 1. 

In Figs. 2 and 3 we present the typical experimental 
time dependencies of L{t) decaying by two-three orders of 
magnitude and analyze them in Sec. IIB. These results, 
partially published in®’^^ are also summarized in Tab.I. 
We demonstrate that the initial stage of decay in all three 
types of the flow, including the coflow regime, typically 
follows a form 

C{t)^ Co = C{t = Q), (3a) 

t-Ti 

with two fitting parameters: the virtual origin time 
Ti and the dimensionless parameter bi. The asymp¬ 
totics (3a) may be rationalized in the framework of the 
Vinen evolution equation for £(t)^® as the decay of ran¬ 
dom tangle of quantized vortex lines. The energy spec¬ 
trum of this vortex tangle in the superfluid component 
(the energy distribution between scales, presented in the 
fc-space) S^{k) is dominated by the intervortex scales 
I. 

The late stage of the decay, discussed in Sec. IIB, fol¬ 
lows a t“®/^-asymptotics^® and may be described by: 

C{t) ^ C^\t) = , (3b) 

[t - 

with two new fitting parameters: T 2 and h^. It is com¬ 
monly believed that the dependence Eq. (3b) is caused 
by the classical Richardson-Kolmogorov cascade in the 
superfluid component, with £g^{k) oc k~^/^ spectrum. 

A natural way to rationalize these observations is to 
assume that the resulting form of the turbulent energy 
spectrum of the superfluid component Eflk) consists of 
two parts: 

£flk) = £’^{k)+£l\k), (4) 

(i) the classical region spanning large scales from the in¬ 
tegral length-scale A down to the intervortex distance £, 
where it is followed by (ii) a quantum contribution (fc), 
corresponding to the random tangle of quantized vortex 
lines, having a form of a peak. Qualitatively, this energy 
spectrum is sketched in Fig. 4(a). 

Section III is devoted to the theoretical discussion on 
the stationary energy spectra £s{k) in the superfluid com¬ 
ponent of “^He quantum turbulence. The main aim of 
this Section is to clarify the striking difference between 
£s{k) in the coflow from that in the counterflow and su¬ 
perflow. In Sec. Ill A, we explain why the energy spec¬ 
trum of channel turbulence in coflow may include, besides 
the classical Kolmogorov spectrum £ (fc) oc the 

quantum peak centered at wave-numbers about 1 j 1. Ad¬ 
ditionally, we argue that this peak has large k asymp¬ 
totics oc I/fc, originating from the superflow near the 
cores of quantized vortices, while for smaller k it bor¬ 
ders the classical spectrum oc k^, caused by the ther¬ 
modynamical equilibrium with the energy equipartition 
between degrees of freedom, see Fig. 4(a). In Sec. IllB 
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we argue that for stationary regimes of the counterflow 
and pure superflow, the energy spectrum of the super¬ 
fluid component has a qualitatively different form from 
that in the coflow, as sketched in Fig. 4(b). The physi¬ 
cal reason for that is the decoupling of the normal and 
superfluid turbulent velocity fluctuations, caused by the 
non-zero value of the counterflow velocity U^s (1) in the 
counterflow and superflow, while in the coflow Uns = 0. 
This position is strongly supported by recent analytical 
researchreviewed in Sec. IllBl. A short discussion 
of a possible form of the energy spectrum, also following 
Ref. is given in Sec. IIIB 2. 

At this point, in Sec IV we are armed to formulate, 
to discuss and to compare with experiments two models 
of decaying quantum turbulence. First of all, based on 
explained above forms of the energy spectra for coflow. 
Fig. 4(a), we suggest in Sec. IV A a “basic model” of 
the VLD decay, Eq. (25), in which the quantum decay 
term (2 la) and the classical energy source term(24e) are 
present from the very beginning of the decay, t = 0. This 
model reproduces both quantum and classical asymp¬ 
totics, given by Eqs. (3), in agreement with the obser¬ 
vations in coflow, shown in Figs. 2 and 3. Based on the 
basic model, in Sec. IV B we present a more detailed (but 
still preliminary) discussion of the underlying physics, 
hidden in its fitting parameters. Although the model 
explains small- and large-time asymptotics of the C{t) 
dependence, it fails to describe the crossover regime be¬ 
tween them, even for the coflow. 

The most striking disagreement between the simple 
“basic” model and observations, seen in Fig. 3, is a 
“bump” (non-monotonic behavior) in the C{t) depen¬ 
dence, in the counterflow and superflow cases. This be¬ 
havior may be naturally explained by the delay in the de¬ 
livery of the energy flux from the classical to the quantum 
part of the spectra (as it has been shown in a slightly dif¬ 
ferent manner in the recent Letter by Gao et for the 
special case of decaying thermal counterflow) required for 
evolution from the more localized in k spectrum, sketched 
in Fig. 4(b), toward the K41-spectrum, shown in Fig. 4(a) 
and (c). This delay is analyzed numerically in Sec. IV C 
and is accounted for in the “improved model” of the 
VLD decay, formulated in Sec. IV D. In this Section we 
show that /l(t)-dependence, following from the improved 
model, allows us to rationalize the main experimental ob¬ 
servations including the bump in the £(t)-dependence at 
the crossover times. 

Following analysis of our experimental findings, we 
suggest in Sec. V an analytical theory of the energy spec¬ 
tra for the steady-state quantum turbulence in the coun¬ 
terflow and pure superflow, based on the algebraic ap¬ 
proximation for the energy fluxes over scales. The re¬ 
sulting spectra are strongly suppressed by the mutual 
friction, leading to the energy dissipation at all scales, 
enhanced by the counterflow-induced decoupling of the 
normal- and superfluid velocity fluctuations. 

In the final Sec. VI we summarize our main results and 
discuss issues that remain out of the scope of this re¬ 


search. In particular, we stress that our simple analyti¬ 
cal theory of the steady-state energy spectra of quantum 
turbulence in the counterflow adopts some uncontrolled 
approximations and simplifications, widely used in the 
studies of classical hydrodynamic turbulence. In the fur¬ 
ther studies of quantum turbulence these assumptions 
have to be either better justified or relaxed. Neverthe¬ 
less, we consider our findings as a natural and perhaps 
even required step in this direction. 


II. DECAY OF THE VORTEX TANGLE IN 
COFLOW, COUNTERFLOW AND SUPERFLOW 


A. Experiment 


The superflow and the coflow of superfluid "^He are 
both mechanically forced by a low temperature bellows 
through a square cross-section brass channel, illustrated 
in Fig. 1. The same channel can be configured to attain 
superflow, coflow and thermal counterflow. For super- 
flow, sintered silver filters - superleaks - block the vis¬ 
cous normal component. Two vertical brass flow chan¬ 
nels have been used, both with the test section 105 mm 
in length and has an internal square cross-section of side 
7 mm and 10 mm, therefore with a factor 2 change in 
cross-sectional area. For coflow, the superleaks are re¬ 
moved and the lower one is replaced by a flow condi¬ 
tioner made from a dense pack of 10 mm long capillaries 
of 1 mm diameter, intended to cut larger scale turbulent 
eddies. The counterflow is prepared by closing one end 
of the channel with a resistive wire heater and leaving 
the other end open. A full description of the mechan¬ 
ically driven superflow apparatus and the measurement 
technique is given in Ref.^"^. Counterflow is studied as in 
previous Prague experiments^®’^®. 

Turbulence is detected by measuring the extra attenu¬ 
ation of second sound caused by the scattering of normal- 
component thermal excitations by the vortex lines. Sec¬ 
ond sound is generated and detected by a pair of vibrat¬ 
ing porous membranes located in the walls of the chan¬ 
nel at its mid-point; the second sound travels across the 
channel, which acts as a resonator. The time dependent 
attenuated amplitude of second sound at resonance a{t) 
can be related to the instantaneous VLD £{t) (assuming 
random and not extremely dense tangle^^) through the 
equation: 


m 


SttA/o r oq 
Bk l-a(t) 


( 5 ) 


where oq and A/o are the amplitude and full width at 
half maximum of the second sound amplitude resonant 
curve for quiescent helium, and B is the mutual friction 
coefficient of order unity, tabulated in Ref.^^. The at¬ 
tenuation of second sound measures the length of vortex 
line per unit volume weighted by a factor sin^ 0, where 9 
is the angle between any element of vortex line and the 
direction of propagation of the second sound. 
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2 
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4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

# 

type of 

T, 

Ua, 

Fn 

t/ns 

£0 • 10 -^ 

Rer 

61 

Tl 

b2 

T 2 

£o/£i 

C 1 /C 2 

tl 

t2 

di 

do 

Col Co 


the flow 

K 

cm/s 

cm/s 

cm/s 

cm“^ 

- 

- 

s 

- 

S 

- 

- 

s 

S 

- 

- 

- 

1 



0.69 

0.69 

0 

0.86 

139 

-1.27 

-1.26 

- 2.22 

-0.08 

2.3 

3.9 

2.0 

14 

0.54 

-0.0035 

8.97 

2 


1.35 

4.99 

4.99 

0 

13.3 

797 

-1.19 

-0.26 

- 0.22 

-0.05 

5.9 

4.3 

1.4 

5.3 

1.73 

-0.019 

25.4 

3 

co¬ 


22.4 

22.4 

0 

106 

3092 

-1.15 

- 0.11 

- 0.001 

1.15 

12.5 

4.7 

1.3 

3.6 

5.36 

0.0049 

58.8 

4 

flow. 


- 

- 

0 

6.0 

843 

-1.04 

- 0.8 

-0.03 

- 1.0 

- 

- 

- 

1.7 

3.05 

-0.0092 

4.16 

5 

7mm 

1.45 

4.99 

4.99 

0 

20.0 

3276 

-1.09 

-0.28 

- 0.05 

- 0.2 

5.6 

3.6 

1.6 

4.5 

3.70 

-0.0039 

20.2 

6 



22.4 

22.4 

0 

80.0 

- 

- 1.25 

-0.13 

0.006 

0.5 

7.7 

13 

1.2 

5.0 

7.92 

0.0046 

100 

7 

counter- 


-0.09 

0.95 

1.05 

0.96 

0.96 

- 1.02 

-0.41 

0.03 

6.0 

2.0 

7.0 

1.5 

2.7 

0.24 

0.0547 

14.0 

8 

flow. 

1.45 

-0.29 

2.89 

3.18 

11.2 

775 

- 0.90 

-0.04 

0.18 

0.1 

6.3 

2.9 

0.7 

8.1 

0.25 

0.0070 

18.3 

9 

10 mm 


-0.77 

7.77 

8.54 

92.5 

1888 

-0.61 

-0.013 

-0.25 

- 0.01 

13.2 

38 

0.1 

3.1 

0.45 

-0.0019 

502 

10 

super¬ 


0.96 

0 

0.96 

1.21 

176 

- 1.40 

-0.18 

- 0.01 

-7 

6.3 

2.9 

0.4 

32 

0.18 

-0.0340 

18.3 

11 

flow. 

1.45 

2.67 

0 

2,67 

10.4 

467 

0.49 

0.07 

-0.003 

-5.8 

17.2 

3.9 

0.2 

6.0 

0.22 

-0.0516 

67.1 

12 

10 mm 


7.41 

0 

7.41 

113 

1246 

0.25 

0.04 

-0.015 

-0.15 

9.1 

5.5 

0.1 

12 

0.68 

-0.0065 

50.1 


TABLE I: Column ^2 shows types of the flow and channel width of chosen twelve sets of experiments, numbered in first column 
from 1 to 12). Next 17 columns display: #3 - Temperature T in K; ^4, ^5 and #6 - superfluid, normal fluid and counterflow 
velocities Ua, Un and Uns in cm/s; #7 - initial VLD Co in cm“^; ^8 - Reynolds number Rer, estimated via counterflow velocity 
Una and normal-fluid kinematic viscosity Vn’, ##8-12: parameters of the fits (3a) for the Vinen’s ( oc t~^) decay and for the 
hydrodynamic (oc decay(3b), namely #9 and #11 -- dimensionless &i and 62 and #10 and # 12 - time-origins ti and 

T 2 in seconds; #13 - ratio of the initial VLD Co to Ci, lowest value of C{t) in the t“^-fit at the time t = t\ (shown in #15); 
#14 - ratio of the VLD Ci to T 2 , the initial (largest) value of C{t) in the t“®^^-fit at the time t — f 2 (shown in #16); #17 and 
#18 - model parameters di and 6 , 2 , related to 61 and 62 by Eqs. (22) and (26b) (with a = 0.06 for T = 1.45 K and a = 0.04 for 
T = 1.35 K). #19 - ratio of the initial VLD Co to C 2 - 


We have performed mechanically-driven coflow, su¬ 
perflow and thermally driven counterflow decay mea¬ 
surements, in the two channels of widths 2A = 7 mm 
and 10 mm, in the temperature range between 1.25 K 
to 2.10 K, and for different values of velocities chosen 
such to produce initial steady-state VLD Cq, spaced al¬ 
most exactly one decade apart: 10"^ cm“^, 10^ cm“^ and 
10® cm“^. For every combination of temperature and 
starting line density we have measured typically 150 in¬ 
dividual decays, under nominally identical experimental 
conditions and we have ensemble-averaged these samples, 
allowing us to resolve up to 4 orders of magnitude of de¬ 
cay on C{t). 


B. Experimental data and their preliminary 
analysis 

In this paper, we restrict ourselves by discussing 12 
typical data sets of experiments, characterized in Tab. I, 
which are chosen to illustrate the basic ideas of this 
research. For T = 1.45 K we discuss three sets with 
£0 ~ 10^ ) 10^ and 10® cm“^ for coflow (sets # 3, 4 and 
5), counterflow (sets # 7, 8 and 9) and superflow (sets # 
10, 11 and 12). In addition, we analyze the coflow sets 
for lower T = 1.35 K with the same Lq ~ 10^ , 10® and 
10® cm“^ (sets # 1, 2 and 3). 

The experiments were conducted in the channels of 
two different widths. For coflow we have more represen¬ 
tative data taken from experiments with 7 mm channel. 


while for superflow and counterflow more representative 
data were obtained with 10 mm channel. We found no 
significant differences between 7 mm and 10 mm channel 
results for a given flow type and therefore compare below 
the most representative results regardless of the channel 
width. 


1. Coflow 

We begin to analyze the time evolution oiC{t) from the 
decaying coflow turbulence, which, in some sense, is the 
simplest case. The VLD C{t) is decaying monotonically, 
as illustrated in Fig. 2 for both temperatures T = 1.35 K 
[shown in Figs. 2(a) and 2(b)] and T = 1.45 K [shown 
in Figs. 2(c) and 2(d)] and for all three initial values of 
£ 0 : 10^ cm“^ (blue lines), 10® cm“^ (green lines) and 
10® cm“^ (red lines). 

There is no qualitative difference between two pre¬ 
sented temperatures, except that for higher T = 1.45 K 
the VLD £(t) is decaying slightly faster. Thus there is 
no reason to analyze these cases separately. 

Our first step is an analysis of the initial stage of the 
coflow turbulence decay shown in Fig. 2(b) and Fig. 2(d) 
(for T = 1.35 K and T = 1.45 K) fitted by t“^-law (3a), as 
shown by dashed lines. The fitting parameters bi and ti , 
given in Tab. I, will be discussed later. For now, we notice 
that the negative virtual origin time ri increases (and 
becomes closer to zero) with increasing £ 0 , as expected 
from Eq. (3a). 
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(a) (b) 






FIG. 2: (color online) Experimental data of the VLD decay C{t)/C q in the coflow (in the 7 mm channel) , normalized by initial 
VLD Cq. The lines correspond (from top to bottom) for Co ~ 10"^ (blue lines), 10® (green lines) and 10® (red lines). The 
explicit values of Co are shown in figures. Data for T = 1.35 K are shown in panels (a) and (b), and for T = 1.45 K - in panels 
(c) and (d). The panels (b) and (d) emphasise the details of the short-time behavior. Quantum t“^-fits, Eq. (3a), are shown 
by dashed dark green lines, while the classical t“®^^-fits, Eq. (3b) - by black dot-dashed lines. 


As seen in Figs. 2(b) and 2(d), the “quantum” 
fit (3a) agrees with the experimental data over the time 
interval 0 < t < about 1 4- 2s, when C{t) decays from 
£o to C{ti) = C\. Ratios of the initial and final values 
of the VLD in the quantum decay, C^jCi and final time 
ti, are given in Tab. I. One sees that for the largest £o ~ 
10 ® cm“^ the ratio £o/£i reaches one order of magnitude. 
On the other hand, sometimes (e.g. set # 4 with T = 
1.45 K and smallest £o — 6 'the quantum regime of 
the coflow decay does not show up. The consequences of 
these important observations will be discussed in Sec. IV. 

The second step is the analysis of the later stage of 
the coflow turbulence decay shown by dash-dotted lines 
for T = 1.35 K and T = 1.45 K in Figs. 2(a), 2(c) and 
fitted from time t 2 by t“®/^-law (3b). In Tab. I we present 


the fitting parameters 62 and T 2 , the starting fit-time ^2 
together with the ratio of the final value £1 = £(ii) in 
the quantum fit to the initial value £2 = £(^ 2 ) in the 
classical fit. 

The intermediate regime between the quantum and the 
classical ones (lasting from ti to t 2 and during which 
£(t) decays from £1 to £ 2 ) is the subject of a separate 
discussion in the Sec. IV. 

Returning to the classical stage of the decay, notice 
that similar to the case of quantum decay, the origin time 
T 2 for £0 ^ 10 ® is larger (closer to zero) than that for 
smaller £0 ^ 10"^, as expected from Eq. (3b). However, 
for larger £0 ^ 10 ®, the time T 2 becomes positive, which 
formally contradicts Eq. (3b). This can be explained as¬ 
suming that there are some transient processes between 
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(a) 


(b) 






FIG. 3: (color online) Experimental data for the VLD decay C{t)/C q at T = 1.45K, normalized by initial VLD Cq. The results 
for counterflow are shown in panels(a) and (b), for superflow - in panels(c) and (d). In all panels the lines correspond to (from 
top to bottom) Co ~ 10"^ (blue lines), 10® (green lines) and 10® (red lines). The explicit values of Co are shown in the figure. 
As in Fig. 2, quantum t“^-fits, Eq. (3a), are shown by dashed dark green lines, while the classical t“®^^-fits, Eq. (3b) - by black 
dot-dashed lines. The panels (b) and (d) show the initial stages of the decay. 


quantum and classical stages of the decay that become 
significant for large initial values of £o • A possible phys¬ 
ical reason for this will be suggested in Sec. IV, after dis¬ 
cussion of the steady-state energy spectra of quantum 
turbulence in the coflow. 

Comparing the rates of the decay at early (quantum) 
and late (classical) stages in Figs. 2, one sees that classi¬ 
cal decay is slower than the quantum one. The physical 
reason for that is simple: in the quantum regime the 
energy source in the decaying vortex tangle is relatively 
small with respect of its dissipation rate and may be ne¬ 
glected. Later, when the dissipation rate, proportional to 
becomes smaller, the Richardson-Kolmogorov en¬ 
ergy cascade toward small scales, serving as an additional 
energy source for the vortex tangle, should be taken into 


account. 

From the formal point of view, it looks contradic¬ 
tory that the decay is slower than . This can 

be resolved by accounting for the interplay of the ori¬ 
gin times in Eqs.(3). This means that decay is 

only an intermediate regime, valid up to some large 
time fa, when the increasing in time intervortex distance 
m = i/^/m (X becomes of the order of the largest 
scale A. We cannot observe this final stage of the decay 
with £(t) ~ 1 cm“^, due to a large noise level. 

Nevertheless, and this is important for discussion in 
Sec. IV, the classical (3b) describes the observed 

decays during tens of seconds after < 2 , when £(t) de¬ 
creases by more than two orders of magnitude. This is 
seen, e.g. in Fig. 2(a) for set ^ 3 with T = 1.35 K and 
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£o ~ 10 ® cm 


2. Counterflow and pure superflow 

Typical examples of decaying VLD in superfluid 
■^He after switching off the counterflow are shown in 
Figs. 3(a,b) and after switching off pure superflow - in 
Figs. 3(c,d). These data, obtained in 10mm channel for 
T = 1.45 K, are very similar to the data from 7 mm 
channel (not discussed in this paper). The lines are 
marked as for the coflow, Figs. 2, according to initial 
VLD : Co ^ 10"^ (blue lines), Cq ~ 10® (green lines) 
and Cq ~ 10 ® (red lines). 

Similar to the coflow, the decay of C{t) the may be 
divided into three stages: (i) The initial quantum stage 
that agrees with the t“^-fit (3a), in the time interval 0 < 
t < ti, during which C{t) monotonically falls from about 
Co to Cl- The htting parameters bi and ti together with 
the ratios CojC\ and time ti are given in Tab. 1. (ii) The 
final classical stage that agrees with the monotonic decay, 
described by t“®/^-fit (3b). It starts at time t ~ t 2 from 
the VLD £2 and lasts for several tens of seconds. Table I 
presents parameters 62 and T 2 together with the ratios 
C\lC 2 and time t 2 - (iii) The intermediate stage between 
the quantum and the classical one, that lasts from ti to 
t 2 and during which C(f) decays from Ci to £ 2 - The 
most striking feature that qualitatively differs this stage 
in counterflow and pure superflow from that in coflow is 
the non-monotonic character of the decay, clearly seen in 
Figs. 3, especially for large £ 0 . 

Comparing Figs. 2 for coflow with Figs. 3 for counter- 
flow and pure superflow, one sees that there are no qual¬ 
itative differences between all these flows at the initial 
quantum and final classical stages of the monotonic de¬ 
cay of VLD. The same conclusion follows from Tab. I, 
which demonstrates only a small quantitative difference 
between these flows at these stages, perhaps with some 
scattered values of displayed parameters for all flows. 

As we discussed above, the counterflow and pure su¬ 
perflow qualitatively differ from the coflow only at the 
intermediate stage, demonstrating non-monotonic decay. 
As we explain below, this difference is a consequence of a 
very different character of the steady-state energy spectra 
(energy distribution between scales) in the counterflow 
and pure superflow with non-zero values of the counter- 
flow velocity C/„s 7^ 0 from the stationary spectrum in the 
coflow, for which C/ns = 0. The energy spectra of quan¬ 
tum ^He turbulence for all three types of flows, discussed 
in our paper, are the subject of Sec. III. 


III. OVERVIEW OF THE ENERGY SPECTRA 
OF QUANTUM ^He TURBULENCE 

In order to rationalize the discussed above time evo¬ 
lution of the VLD £(C) in counterflow, coflow and pure 
superflow turbulence we need to clarify the energy spec¬ 


tra in all these types of quantum turbulence. We be¬ 
gin in Sec. Ill A, with discussion of the energy spectrum 
of coflow grid turbulence. After that, following recent 
Ref. we overview in Sec. Ill B the effect of the coun¬ 
terflow velocity on the steady-state energy spectra of the 
counterflow and pure superflow turbulence, important for 
the understanding of the appearance of a bump on the 
time dependence of C{t) at the intermediate stage of the 
decay. 


A. Steady-state coflow grid turbulence 

In coflow turbulence, it is very natural to assume that 
the mean velocity profiles of the normal and superfluid 
velocity, Un{y) and Us{y), practically coincide almost ev¬ 
erywhere in the channel, perhaps except for narrow re¬ 
gions near the walls. The reasons are the following: The 
pressure drop is the same for both fluid components, rela¬ 
tively large mutual friction tries to lock the mean normal 
and superfluid velocities and only kinematic viscosities, 
important in the near-wall region (viscous and buffer lay¬ 
ers) are different: t'n 7 ^ t's = 0. This causes some decou¬ 
pling of the velocities in that region. As a consequence, 
the kinetic energy is transferred from the mean flow to 
the turbulent velocity fluctuations. Its value is propor¬ 
tional to the shear of the mean velocity S{y) = dU{y)/dy 
and is practically the same in the entire channel for the 
normal and superfluid components. Thus we assume that 
the energy spectra of the normal and superfluid compo¬ 
nents practically coincide in the entire energy containing 
and inertial intervals of scales, 

£n(fc) = £s(fc). ( 6 a) 

Ignoring unessential for present discussion of intermit- 
tency effects, we can use classical K4I energy spectra for 
both "^He components: 

S:\k) = Sfik) cs , ( 6 b) 

where £„ = £3 = £> are the corresponding energy fluxes 
over inertial interval of scales and C'k 4 i ~ 1.4 is the Kol¬ 
mogorov constant. 

In the classical channel flow the mean turbulent ki¬ 
netic energy and the kinetic energy of the mean flow 
= {U)^ /2 are related by: 

4 

~_ - _ f7) 

[«K'ln(Re,) + Rci ]2 ■ 

Here Reynolds number Re,- = A^/p'A/r'„, where A is 
half channel width, p' is the pressure gradient that drives 
the flow, Kk ~ 0.41 is the von Karman constant and the 
constant Bc\ « 5. Using this equation, we can find 
which is dominated by the outer scale of turbulence ~ A. 
Then, using Kolmogorov energy spectrum ( 6 b), we esti¬ 
mate its part, Ei, originating from the intervortex scale 



t. El ~ . This part corresponds to the en¬ 

ergy of the random vortex tangle of VLD £, that may 
be found from the relation Ei ~ By equating these 
expressions we can roughly estimate VLD connected 
to the classical energy spectrum for different values of 
Us- The result is the following: for smaller [ 4 , that cor¬ 
responds to £ ~ 10^cm“^, the estimated value of VLD 
£* is close or above the experimental value of initial VLD 
Cq: £» ~ (0.9-^ 1.4) £ 0 , while for larger Us (correspond¬ 
ing to £ ~ 10® cm“^) we find that £* ~ (0.5 1.0) £ 0 . 

This means that, in order to explain how the observed 
initial value of VLD £0 may be larger than the classi¬ 
cally generated value £*, one should find some additional 
mechanism of vortex generation, besides the classical flow 
instabilities in the channel flow. 

It is almost common belief (see, e.g.. Ref. ^’^®) that 
large scale motions of normal and superfluid components 
are dynamically locked, in the sense that their turbulent 
velocities coincide everywhere in space at any time: 

u„(r,t) = Us(r,t), (8) 

and thus the counterflow mechanism of the vortex lines 
generation seems to be absent. 

However, as was shown in detail in Refs. the 

normal and superfluid energy spectra deviate from each 
other in the crossover region M ~ 1. The relation (8) is 
also violated in the narrow region near the channel walls 
and near the surface of the grid at the channel entrance, 
where normal and superfluid components satisfy different 
boundary conditions. In these areas the velocities of com¬ 
ponents differ and excite a random vortex tangle leading 
to a peak in the energy spectrum near the crossover scale 
I. 

We thus conclude that the resulting form of the tur¬ 
bulent energy spectrum E(k) may include, besides the 
classical Kolmogorov region (6b) (spanning scales from 
the integral length-scale A down to the intervortex dis¬ 
tance £), also the quantum contribution (k), which has 
a form of a peak corresponding to the random tangle of 
quantized vortex lines. This peak has large-fc asymptote 
oc I/fc originating from the velocity field near the quan¬ 
tized vortex lines. In the region of fc < 7r/£ the quantum 
peak is adjoined by the classical region with the thermo¬ 
dynamic equilibrium spectrum oc describing equipar- 
tition of energy between degrees of freedom®^. 


B. Steady-state counterflow and pure superflow 
^He turbulence 


1. Counterflow induced decoupling in quantum turbulence 

The first important ingredient of understanding of the 
steady-state energy spectrum in counterflow turbulence 
is decoupling of the normal and superfluid velocity fluctu¬ 
ations, Mn(^, t) and Us{r, t), caused by the their sweeping 
in opposite direction by the corresponding mean veloc¬ 
ities and Us- Considering the velocity fluctuations 
of characteristic scale R (referred for the shortness to as 
i?-eddies) one can easily estimate their overlapping time 
as 

Tov {R) = R/Uns , Uns = Un - Us - (9) 

Another characteristic time in the problem of velocity 
coupling is the interaction time Ti„t of normal and su¬ 
perfluid eddies, due to the mutual-friction force acting 
between them. To estimate Tint, we recall®® that normal 
eddies act upon the superfluid eddies with the force (per 
unit mass) Fs 

Fs{r,t) ~ n[Uns + Unir,t) - Us{r,t)], (lOa) 
ft = a{T)K£, 

where is mutual friction frequency for the superfluid 
and a{T) is the dimensionless mutual friction parameter, 
tabulated in Ref. On the other hand, the superfluid 
eddies act on the normal fluid ones with the force 

Fn{r,t) ~ -D„ [Uns + Unir,t) - , (lOb) 

where = flps/pn- As required by the conser¬ 
vation of the mechanical momentum, the weighted 
sum psFs{r,t) + pnFn{r,t) = 0, while the difference 
Fnsirfl) = Fs{r,t) — Fn{r,t) plays an important role 
in the intrinsic dynamics of the velocity differences. In 
the equation for Fns{r,t) 

-^ns(^5 ^) — ^ns[£^ns T ^n(^7 0 f)] 5 ( 10 ^) 

0ns — £^p/Pn — Ons ^£ ^ O^ns — ^(^)p/Pn(£) .(lOd) 

The resulting mutual friction parameter q;„s is only 
weakly dependent on T, varying between 0.7 and 0.5 in 
the temperature range (1.4-^ 1.9) The first term 

in Eq. (10c), iJnsUns, is compensated in the stationary 
counterflow by the force, originating from the temper¬ 
ature drop. The second term, ^}ns[un{r,t) — Us{r,t)], 
responsible for the coupling of the velocities by mutual 
friction, is proportional to the scale-independent “inter¬ 
action frequency” Clns- Accordingly, the mutual-friction 
interaction time can be estimated as follows: 


Detailed analysis of the steady-state energy spectra in 
the counterflow turbulence, that reflects our current un¬ 
derstanding of this problem, will be done below in Sec. V. 
Here we overview only some earlier results on the coun¬ 
terflow energy spectra, which, nevertheless are sufficient 
to formulate in Sec. IV our model of decaying quantum 
turbulence. 


Tint — 




, kC 


( 11 ) 


As was shown in Ref. ^® (see also^®), the overlapping 
time Tov (9) should be compared with the i?-independent 
interaction time (II). If, for large Uns or small R, the 
overlapping time is much smaller than the interaction 
time, Tov ^ Tint, the mutual friction does not have 
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(a) Coflow 


(b) Counter- and super-flow 


(c) Late time asymptotics 





FIG. 4: (color online) Sketch of the stationary superfluid turbulent energy spectrum in log-log coordinates, \og£s{k) vs. \og{k£). 
According to Eq. (4) spectrum £s{k) consists of classical and quantum £^{k) parts, colored in light blue and cyan. In 

coflow (Panel a) £s^(k) consists of cascade part (fe) oc (for k < kx ) and thermodynamic equilibrium part (k) oc k^ 

for 1/1 ^ k > kx- In counterflow and pure superflow (Panel b) the quantum contribution £^ (k) and the classical thermal 

TD 

bath part £^ (fe) look similarly to that in coflow, while the cascade part, supercritical LNV-spectrum, £1 lnv given by Eqs. (19), 
ends at some fe* < l/£ and does not provide energy to the quantum vortex tangle in stationary regime. After switching off the 
counterflow or pure superflow, the spectrum shown on Panel b evolves to that shown in Panel c, switching on the energy flux 
toward quantum vortex tangle after some delay. 


enough time to couple the normal and superfluid velocity 
fluctuations. In this case one expects that the velocities 
remain uncoupled. Otherwise, i.e., for large scales R or 
small C/ns, Ajv 'Lnt, and one should expect full cou¬ 
pling (8). 

Therefore the coupling-decoupling process is governed 
by the dimensionless “decoupling” parameter 

(^(fc) = ~ ^ ^ ^ tt/R. (12a) 

7”ov i ^ns 

The analytical theory of the coupling-decoupling 
processes^^, using Langevin inspired approach to model 
the nonlinear term, results in the analytical expression for 
the cross-correlation Snsik) in terms of £sik) and £-n(k), 
defined by Eqs. (B5). Using Eqs. (13) and (16) from 
Ref. the result in the relevant for our experiments 
range of parameters can be written as: 


fns(fc) 


D{k) 


Ufn(fc) + ^Yv£s{k) 

-u- 

^ns 

^arctan[C(fc)] . 




n + u„,(i2b) 
(12c) 


In Eqs. (12) the dimensionless “decoupling function” 
D{k) is defined via the dimensionless “decoupling param¬ 
eter” C,{k), and describes the decoupling of the normal 
and superfluid velocity fluctuations, caused by the coun¬ 
terflow velocity. The plot of U as a function of k/kx is 
shown in Fig. 5, where 


kx — 2 Uns/Uns, (12d) 


smallest one, ko ~ tt/A, up to kx , which is slightly larger 
than ko and are decoupled in the rest of fc-range up to 
tt // » kx ■ 


2. Qualitative analysis of the counterflow energy spectrum 


The next question of principal importance for the anal¬ 
ysis of our experiments is how the decoupling of the ve¬ 
locities at small scales affects the steady-state energy 
spectra before switching off the flows. We will discuss 
this problem in details later in Sec.V. Nevertheless, to 
get a qualitative form of the energy spectra, sketched 
in Fig. 4, required for formulation of our model of de- 



The analysis of Eqs. (12), given in Ref. allows us to 
conclude that the normal and superfluid velocity fluctu¬ 
ations are coupled in the wave number interval from the 


FIG. 5: (color online) Decoupling function D{k/kx) = 
£ns/£mi VS dimensionless wave number k/kx- 
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caying quantum turbulence, it is sufficient to overview 
relatively resent results on energy spectra in turbulent 
^He, see Refs.^^’^^ and analysis, presented in Ref. 

The coupling - decoupling process of normal and su¬ 
perfluid velocity fields crucially affects the energy dissipa¬ 
tion due to mutual friction and, correspondingly, modifies 
the energy spectra Ss(k,t) and To demonstrate 

this, we use the evolution equations for these objects, 
derived in Ref. : 

= £l[£Uk,t) - £s{k,t)\ , (13a) 
+A/'^n = n^[£nsik,t) - £nik,t)] .(13b) 


Volovik (LNV) in current notations is: 



n k \2/3 
flcr ^ ^kn) 


n i 2 

Her - 


, (16a) 


(16b) 


At n = ricr it becomes the scale-invariant “critical” spec¬ 
trum 


4r(fc)=£o(y)', 


LNV-critical. 


(17) 


For < flcr solution (16a) can be considered as “suberit- 
ical” and written as follows: 


Here AfCg^n are nonlinear terms, which conserve kinetic 
energy and therefore may be presented in the divergent 
form: 


A/-/:. = ^ , = ^ . (13c) 

Here £s{k) and £n(fc) are the energy fluxes in correspond¬ 
ing subsystems. For k ^ fcx, due to the decoupling, 
£nsik) <C £s{k). Therefore the term £nsik) may be ne¬ 
glected in the RHS of Eq. (13a), which becomes —flgSg. 
This is similar to the equation for £s for the quantum 
turbulence in superfluid ^He-B, written below for the 
steady-state case: 


= 0. (14a) 

dk 

Notice that the energy flux over scale £{k) can be ex¬ 
pressed in the terms of the third-order correlation func¬ 
tion of tts(fc) exactly as in classical turbulence. In order 
to proceed further, one can borrow a closure procedure 
from classical turbulence that expresses Ssik) in terms 
of the energy spectrum £sik). Even though this step 
is widely used, it is worth remembering that it is uncon¬ 
trolled. The simplest algebraic closure relation suggested 
by Kovasznay^® 

(14b) 

just follows from the K41 dimensional reasoning 

£si£\k) = , Ck = (S/d)^/^ « 1.4 . (14c) 

The prefactor | is chosen to simplify the appearance of 
some of the equations below and to get numerical value 
of the Kolmogorov constant Ck reasonably close to its 
experimental value. 

Ordinary differential Eqs. (14) 

l^[k^/^£!/^{k)]+n£M = 0, (15) 

should be solved with the boundary conditions £s{ko) = 
£o at the lowest wave number ko in the inertial inter¬ 
val. The physical solution, found by L’vov, Nazarenko, 


1-3 p 1 

4b(fc) = 4 ° ^ 


1 


kcT — ki 


kb/3 Lfc2/3 ^2/3 
/ n \3/2 
V Opr - H / 


, LNV-subcritical, 

(18a) 


If O is close to Opr from below. Opr — O <C Opr, and k 
is smaller than the crossover wave number fcpr, 4b(A) is 
close to the critical solution (17). However at k /S> kcr, 
the subcritical spectrum (16a) has K41 asymptote 


(- 1 ;)' 38 .) 

with the energy flux Soo smaller than the energy influx 
Eq. The difference Eq — Eoo is dissipated by the mutual 
friction. 

In the case O > Opr, Eq. (I6a) can be considered as 
“supercritical” spectrum 


£sp{k) = 4 ^ ^ 


I l2 


0^5/3 Lfc2/3 ,2/3 


fc* = k. 


/ O \3/2 
\ O - Opr / 


LNV-supercritical, 

(19) 


that terminates at some k = kn. However, if O is close 
to Opr from above, O — Opr ^ Opr, such that 4, 

and k < /c*, we see that the supercritical spectrum is 
close to the critical one. Examples of three versions of 
the LNV spectra are plotted in Fig. 8. Note in passing 
that similar conclusions on behavior of ^He-B quantum 
turbulence follow from independent analysis of Vinen^®. 

All this means that, assuming full decoupling of the 
velocity fields, the behavior of the superfluid component 
of "‘He in the counterflow becomes similar to that in ^He- 
B turbulence with the normal fluid component at rest. 
Therefore one can expect that the energy distribution 
between scales for k /$> k^ may be described by the LNV- 
spectrum (16) or similar localized spectrum. 

For k < kx, due to partial velocity correlations, the 
energy dissipation due to mutual friction is much weaker 
than for k > kx, although it cannot be completely ne¬ 
glected as in coflowing ^He (with classical K4I energy 
spectrum). Thus we can expect only moderate suppres¬ 
sion of the energy spectrum as compared to the K4I case. 
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as was recently observed in Ref. and analysed below 
in Sec. V. The resulting energy spectrum is sketched in 
Fig. 4(b). Besides the classical part, this spectrum in¬ 
cludes a quantum peak, presumably more intensive than 
the peak in the coflow, shown in Fig. 4(a). These peaks 
have similar structure: large fc-asymptotics ex Xjk and 
small fc-asymptotics oc . 

IV. BASIC AND IMPROVED MODELS OF 
DECAYING QUANTUM TURBULENCE 

In this Section, we first propose and solve a basic model 
of the VLD evolution, L{t), in decaying quantum turbu¬ 
lence. The model will then be further developed to clarify 
the experimental facts in more detail. 


A. Basic model of decaying coflow turbulence 


r;"' ~ e/K^ . (23) 

Integrating Eq. (6b), we And the total classical kinetic 
energy 


■K j Z 

= J E^\k)dk~C^,^e‘^/^ J k-^/^dk 

tt/A 






^£C?2A^2/3 


(24a) 


Here we roughly estimated lower and upper limits of the 
inertial interval as tt/A and 7r/f. Bearing in mind that 
at our experimental conditions A significantly exceeds f, 
we neglected f with respect to A. Finally, we replaced 
in Eq. (24a) all numerical factors by a new fitting dimen¬ 
sionless constant c ?2 — 1- As a result 


Assume that the time derivative of C{t) consists of a 
simple sum of the quantum decay term —rj^n and the 
classical source term rjci, neglecting possible processes of 
their interaction: 


d£{t) 

dt 


_jyQ _|_ jjC 


( 20 ) 


£~2(£:"i)3/2/((i2A) . (24b) 


Now we can write the classical energy balance equation 


dt ^ ^2^ 


(24c) 


The quantum tangle decay (without counterflow veloc¬ 
ity) is usually discussed in the framework of the Vi- 
nen equation with the (quantum) decay term = 
X 2 kC '^Estimating the phenomenological coeffi¬ 
cient X 2 in the local induction approximation^®’^®, we 
rewrite in terms of the dimensionless mutual friction 
parameter a and, the vortex line curvature S, normalized 
by the intervortex distance t', C 2 ^ 1: 

rf^ a K , C 2 = S£, A = ln(f/ao) . (21a) 


with the solution: 


5"' 



d2A 

t - T2 


2 



(24d) 


where £q^ is the initial energy at t = 0. Together with 
Eqs. (24b) and (23) this finally gives 




£^ 


2{d2Ay 

{t-T2Y 


(24e) 


Bearing in mind that for our conditions the parameter 
A/(47r) ~ 1 varies from 0.9 to 1.1 and C 2 — 1 depends 
weakly on temperature^^, we can simplify this relation 
by introducing a dimensionless fitting parameter di , and 
write 


Collecting Eqs. (20), (21b), (24d) and (24e), we finally 
suggest a simple model of the VLD decay in the form: 


dC 2 ((i 2 A)® an 
dt K^(t — T2Y 


TyQ = ^£2 . (21b) 

di 

With this vp and without classical source term 
Eq. (20) has well known U^-solution (3a) with 

bi= di/{a kCoTx) . ( 22 ) 

The positive contribution rf^ in the RHS of Eq. (20) 
originates from the direct energy flux £ from the classical 
energy containing scale kA ~ 1, to the quantum energy 
peak £^{k), located at kd ~ 1, see Fig. 4. To estimate 
ry'^*, recall that ~ (see also Ref. ^®). Therefore 
the flux of VLD from the classical scales to the quantum 
energy peak may be estimated as 


Here di and d 2 are dimensionless phenomenological pa¬ 
rameters. Note that T 2 is of the order of turnover time 
TA — A/^/^ of the largest (energy containing) eddies 
in the flow. 

When the classical energy £g* is small, the virtual ori¬ 
gin time T 2 , according to Eq. (24d), is large and the first 
term in the RHS of Eq. (25) may be neglected. In this 
case Eq. (25) may be easily solved, giving U^-decay 
Eq. (3a). Assume now that for large times the time 
derivative d£/dt in Eq. (25) may be neglected. Then 
Eq. (25) reproduces t“®/®-decay Eq. (3b) with 
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(a) n 


(b) bi 


(c) di 



FIG. 6: (color online) Virtual origin times ri and T2 for the quantum (3a) [panel (a)] and classical (3b) decay [panel (d)[ vs 
initial VLD jCq in initial coflow, counterflow and pure superflow. Panels (b) and (e) display /!o-dependence of the bi and 62 
fit parameters. Panels (c) and (f) display Ho-dependence of the di and d 2 fit parameters. The numerical values of all these 
parameters are given in Tab. I. 


or by using di from Eq. (22) 


d2 = &2 


3/2 


A 



(26b) 


Neglecting for simplicity the time dependence of the 
first term in the RHS of Eq. (25), one can solve it exactly 
with the result: 


C{t) = C^^{t) coth 


-cW)- ■ 


Remind that 


(27) 


lim coth[a;] —> 1/x, lim coth[a:] —)■ 1. 

x—vO x—^00 


Bearing in mind the approximate character of the evo¬ 
lution model Eq. (25), we will consider Eq. (28) for C{t) 
as an analytical form of the VLD time dependence C{t), 
practically equivalent to suggested basic model of decay¬ 
ing superfluid turbulence. 


B. Analysis of the decay fitting parameters 

The £o"dependences of the fitting parameters ri and 
T2, bi and 62, di and ^2 for three different types of the flow 
are given in Tab. I. For clarity, we additionally display 
these dependencies in Figs. 6. 


Therefore,when ^ £{t) —)■ . Otherwise, in 

the limit < C^\ C{t) For ~ the func¬ 

tion C{t) describes a smooth transition between and 
being always larger than both and In this 
way, the approximate solution (27) for L{t), interpolates 
the solution of the model Eq. (25) between the small- and 
large-time asymptotes. 

Even better interpolation between exact asymptotes 
and gives the following modification of Eq. (27): 


C{t) = /l‘^'(t) coth 


r/:‘=^(0)i 


(28) 


1. Quantum and classical origin times ri and T 2 

The dependence of the quantum origin times ti on the 
initial VLD Cq for coflow, counterflow and pure super- 
flow is shown in Fig. 6(a). As expected from Eqs. (3a) 
and (22), the values of ti for coflow are negative and 
monotonically increase with Cq. There values of ti for 
coflow are clearly different from those in the counterflow 
and pure superflow, especially for smaller Cq. This is 
related to the fact that in the coflow, the classical en¬ 
ergy flux towards the quantum vortex tangle is present 
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from the very beginning of the decay (at t = 0), while 
in the counterflow and pure superflow it appears only af¬ 
ter some delay, required for developing the Kolmogorov 
cascade from the initial spectrum sketched in Fig. 4(b). 

As expected, the classical origin time T2 tends to de¬ 
crease for larger Lq, as depicted in Fig. 6(d). The most 
striking fact, seen in Fig. 6(d), is that T2 can be posi¬ 
tive. This definitely contradicts the basic model of £{t)- 
decay, formulated in Sec. IV A. Non-monotonic evolution 
of C{t), clearly seen in Figs. 3, also calls for improve¬ 
ments of the basic model of £(t)-decay, in order to ac¬ 
count for the time delay of the energy flux into the quan¬ 
tum vortex tangle. This is done in Sec. IV D, prefaced by 
Sec. IV C, devoted to the study of the energy-flux delay in 
the decaying turbulence in counterflow and pure super- 
flow, caused by the later development of the Richardson- 
Kolmogorov cascade from the localized energy spectra. 


2. Quantum and classical parameters bi, di and &2, d 2 

The £o-dependencies of the Ht parameters &i, di and 
62, d 2 are shown in Figs. 6(b, c) and Figs. 6(e, f). 

Quantum parameters &i, di are of the order of unity; 
their deviation from unity reflects the non-universal char¬ 
acter of the transient regime after switching off the flow: 
the mean velocity has to relax to zero, anisotropic statis¬ 
tics of the vortex tangle, affected by the mean flow, has 
to become isotropic, etc. We will not discuss these com¬ 
plicated issues in this paper, as they do not seem to con¬ 
tribute to better understanding of the basic physical pic¬ 
ture of decaying quantum turbulence in three types of 
the flow. 

The same can be said about the scatter of values of the 
classical fitting parameters &2, ^2 for different types of the 
flow. Notice only that their values, being much smaller 
than unity, become even smaller for larger Cq. The same 
tendency is demonstrated by the ratio C 2 (from which 
the classical decay begins) to the initial value of VLB Cq. 
This may be easily interpreted in the following way: the 
random vortex tangle decays according to the quantum 
law as long as the energy influx from the classical part 
of superfluid turbulence may be neglected. This holds as 

K41 

long as the tail of classical energy spectrum is 

smaller than the quantum energy spectrum E^{'k/E) at 
this scale. Therefore the ratio {-k/ 1)/E^'^^ (jr/(.) may 
be roughly estimated as 1 /b 2 1 or as a value between 

the ratios C^jCi and C^jC 2 (remind that £1 is the VLB 
at which the quantum decay terminates). As shown in 
Tab. I, both ratios are much larger than unity and both 
tend to increase for larger £o- 

Based on this analysis, we conclude that the quantum 
peak of energy at the intervortex scale £, as a rule, dom¬ 
inates over the tail of classical energy spectrum at this 
scale, as depicted in Figs. 4. Notice that, generally speak¬ 
ing, this qualitative conclusion may be guessed just from 
the observation of the decay laws, shown in Figs. 2. What 


is added by our analysis is a semi-quantitative estimate 
of the ratio E'^ {tt/£) /[t:/E). 

C. Energy-flux delay in the decaying counterflow 
turbulence 

It is commonly accepted that the Richardson- 
Kolmogorov cascade develops from any localized in the 
fc-space initial state over finite time, of the order of a few 
turnover times of the energy-containing eddies. Never¬ 
theless, the details of the transient regime and how they 
depend on the initial state are poorly understood. In 
order to clarify the law of delay in “switching on” the 
energy flux Sci{t), Eq. (24e), that contributes to the RHS 
of the basic model (25), we adopt in our paper so-called 
Sabra-shell model of turbulence, successfully utilized in 
studies of quantum turbulence, e.g., in Refs. 

The required version of the Sabra model and the nu¬ 
merical procedure are described in Appendix A. Here we 
present only results of the time evolution of the energy 
spectrum Es{km,t) = (|umP) /km in the decaying quan¬ 
tum turbulence from five types of the initial conditions 
Es{k, 0) = Eo{k), shown in Fig. 7: 

ICa - K41 energy spectrum (6b), E^^aik) oc 

ICb - Experimental counterflow spectrum, Eo^b{k) oc 
reported in Ref. 

ICc - Critical LNV-spectrum (17) Eo^dk) oc k~^; 

ICd - Supercritical LNV-spectrum (19) Eg^dik) = 0 for 
k> k^.] 

ICe - Subcritical LNV-spectrum (18) Eo^dk). 

The time dependence of the total energy for five types 
of the initial conditions, (a)-(e), are shown in Fig. 8(a). 
As expected, in all cases the large time asymptotics, 
E{t) oc agrees with Eq. (24d) - see black dashed line. 
What is important for current discussion, is the small 
time behavior for t below few ta, shown in the inset 
of Fig. 8(a). Clearly, for critical and supercritical initial 
conditions, ICc and ICd, the energy E{t) (shown by co¬ 
inciding solid red and dashed light blue lines) does not 
decay up to t « ta. This is the time required for devel¬ 
opment of the Richardson-Kolmogorov cascade, transfer¬ 
ring energy into the dissipative range of scales (large k). 

On the contrary, for the K41 initial condition, ICa, 
the energy decays from the very beginning - see blue 
solid line in the inset of Fig. 8(a). This behavior is also 
expected. For two intermediate initial conditions, ICb 
and ICe, the initial energy is also localized in the region 
of small k, but not so strongly, as in ICc and ICd cases. 
Accordingly, the time dependence of E{t) for ICb and 
ICe cases, shown by green and orange lines, demonstrate 
intermediate behavior. 

All these features are clearly seen in Fig. 8(b), showing 
the time dependence of the rate of energy dissipation 
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e(t). In our simulation we used a very small value of the 
kinematic viscosity, therefore e{t) is actually a measure 
of the energy flux via crossover wave number ^ l/£ that 
serves as the energy flux from classical to quantum scale 
range. Thus, according to Eq. (23), we can say that the 
time dependence £{t) actually gives the time dependence 
of the important ingredient of the model Eq. (20), the 
classical energy source ~ e(t)/K^. 

All plots of in Fig. 8(b) have asymptote for 
large t, in agreement with Eq. (24e). Moreover, for the 
K41 initial conditions, ICa, expected in coflow decay, and 
shown by the blue line, Eq. (24e) provides reasonable fit 
of for all times, as shown by the solid black line 

in Fig. 8. Thus our numerical simulations support the 
basic model (25) of the decay in coflow, formulated in 
Sec. IV A. 

The situation is completely different for other initial 
conditions, expected for the counterflow and pure su¬ 
perflow cases. There is the most striking difference for 
the well localized, critical and supercritical ICc and ICd, 
shown in Fig. 8 by (practically coinciding) red and light 
blue lines. One sees that for small time « 0, then 
it is sharply switching on and after few turnover times 
ta reaches the “basic” behavior (24e). This can be ac¬ 
counted for by “improving” the basic model, introducing 
into Eq. (24e) for the classical source term the time delay 
function Fdei(t): 



FIG. 7: (color online) Five variants of the stationary energy 
spectra, compensated by fe®f®i5a(fc), serving as the ini¬ 

tial conditions for the Sabra-decay Eq. (A2): a) K41 energy 
Eo,a oc Eq. (6b) (solid blue line); b) Experimental 

counterflow spectrum, Eo,b{k) oc reported in®^; (green 

line with squares) c) Critical LNV-spectrum (17) Eo^d{k) oc 
k~^ (red line with diamonds); d) supercritical LNV spectrum 
Eo,d{k) = 0 for k > fc*, Eq. (19); (light blue line with circles) 
e) subcritical LNV spectrum (18) (brown line with triangles). 




cl 


2{d2^f 

{t + Tdel - T-2)3 


^ Tdel(t) 


2{d2^f 

(t-T2)3 ■ 


(29a) 


D. Improved model of VLD decay vs experiment 

1. Improving the basic model by the energy-flux delay 


It is convenient to choose Edei(t) as a square of new func¬ 
tion fn{i) that has a simple form: 


Fde\{t) = flit) , fn{t) = ■ (29b) 

^ Mel 

As required, /„(t) —>■ I for t —>• oo. Generally speak¬ 
ing, /(O) = 0 only for the supercritical case, ICd. For 
all other cases (except for the K41 spectrum, which we 
are not discussing here) /(O) 0, although small. This 

difference is not important for us and for simplicity we 
adopted in Eq. (29b) a simple assumption that /n(0) = 0. 
The delay time Tdei in Eq. (29b) is expected to be about 
the largest eddies turnover time, ta. Indeed, as seen in 
Fig. II, Tdel - 0.4ta. 

Notice that the parameter n in Eq. (29b), responsible 
for the sharpness of the delay function, is different for 
different initial conditions. Figure 11 demonstrates that 
for the weakly localized initial conditions, ICb and ICe, 
a reasonable approximation to the numerical observation 
can be reached with n = 1, while for strongly localized 
initial conditions a good fit corresponds to n = 6. 


In the previous Sec. IV C we improved the classical 
source term accounting for the time-delay by the delay 
function Fdei(t). Substituting the new form (29) of 
in the basic model Eq. (25), we formulate the “improved 
decay model” of quantum turbulence: 

(K ^ 2(d2A)^fdel(t) _ ^^2 /gQ', 

dt K^it + Tdel - T 2 )^ di 


For Tdec = 0 the improved model (30) coincides with the 
basic model (25). For t < Tdei the delay function (29b) 
Fdei(t) < 1 and the energy flux term in Eq. (30), that 
is proportional to E(t), is suppressed. For t ~ Tdei this 
term is gradually switching on and, finally, for t ^ Tdei 
the improved and basic models, Eqs. (30) and (25) again 
coincide. 

By analogy with Eq. (28), we can formulate the ana¬ 
lytical form of the improved decay model: 


Lit) = /dei(t)>C“^(t + Tdel) coth 


r/dei(t)/:“no)i 

AQ 


(31) 


Now we are fully armed to compare the suggested an¬ 
alytical models (28) and (31) with experimental observa¬ 
tions. This comparison is the subject of the two following 
Sections. 
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(a) (b) 




FIG. 8: (color online) Panel (a) - time dependence of the total energy £{t) in decaying superfiuid turbulence from different initial 
conditions, shown in Fig. 7. The inset shows the short-time evolution of the total energy S(t). Panel (b) - time dependence of 
the rate of energy dissipation e{t)/K^ that plays role of the VLD flux in the model Eq. (20). The lines are colored as in 

Fig. 7 without markers. In panel(a), the lines for critical and subcritical LNV-spectra coincide. The lines for experimental and 
supercritical LNV-spectrum almost coincide. 


2. Basic decay model vs coflow experiment 



FIG. 9: (color online) The details of the short-time depen¬ 
dence of the rate of energy dissipation e(t)/K^, shown in 
Fig. 8(b). Black solid line - the basic model dependence (24e) 
for r?‘^*(t) fits K41 initial conditions IGa. Black dashed line 
shows improved model dependence (29) with n = 2; it ap¬ 
proximates numerically found r?‘^*(t) for weakly localized ini¬ 
tial conditions ICb and IGe. Black dot-dashed line, given by 
(29) with n = 6, fits r?'^*(t) computed with strongly localized 
initial conditions ICc and IGd. 


We begin here with the more simple coflow case, choos¬ 
ing for comparison the less noisy set ^ 3, demonstrating 
(red line) in Fig. 2(a) the decay of almost 4 decades over 
50 s. This line is reproduced (in red) in Fig. 10(a) to¬ 
gether with the plot of analytical C{t) (shown by the blue 
dotted line), predicted by the basic model. This line is 
barely seen because it practically coincides with the red 
experimental line within the line-width, broadened after 
25 s by noise. Some discrepancy between the experiment 
and the model prediction is better seen on the closeup 
Fig. 10(b) showing the first 5 s of the decay. 

Notice that the basic model is very simple: it com¬ 
pletely ignores any interaction between the quantum 
peak and the classical large-scale turbulence. In formu¬ 
lating the model we also ignored the energy stored in the 

TD 

equilibrium part of the classical energy . Last but not 
least, the model does not account for the spatial inhomo¬ 
geneity of the turbulent channel flow, in which turbulent 
kinetic energy significantly depends on the distance to 
the walls. Bearing all these in mind, we consider the 
agreement between the experiment and the model quite 
satisfactory. 

3. Improved decay model vs counterflow experiment 


Our next step is a discussion of the decay in the coun¬ 
terflow and pure superflow cases. For our rather sim¬ 
plified manner of modeling, there is no physical differ¬ 
ence between the counterflow and pure superflow decay¬ 
ing turbulence. Therefore from the data presented in 













16 


(a) 


(b) 




FIG. 10: (color online) Comparison of the experimental observation of cofiow decay with the basic model predictions. Red solid 
lines, reproduced from red lines in Figs.2(a,b), show the coflow data for set ^ 3 (T = 1.35K, £o ~ 10®cm“^). Blue dotted 
lines show basic model prediction, £(t ), Eq. (28) with ri = 0.5s, T 2 = Is and 62 = 0.075. 


(a) 


(b) 




FIG. 11: (color online) Comparison of the experimental observation for the superflow with the improved model predictions. 
Red solid lines, reproduced from red lines in Figs. 3(c,d) for superflow # 12 (T = 1.45K, jCq « 10®cm“^). Blue dotted lines 
show improved model prediction, jC(t ), Eq. (31) with n = 2, ti = 0.015 s, T 2 = 0.1s, Tdec = 0.6 s and 62 = 0.22. 


Figs. 3, we chose the set # 12 (superflow at T = 1.45 K, 
£0 — 10®cm“^) shown in Figs. 3(c,d) by red lines. In 
this case the noise is relatively low and the bump on 
the £(t)-dependence is clearly pronounced. Red lines in 
Figs. ll(c,d) reproduce the experimental results shown in 
Figs. 3(c,d) in red. Blue dotted lines result from the an¬ 
alytical form £(t) of the improved model of decay, given 
by Eq. (31) with n = 2, ti = 0.015 s, T 2 = 0.1 s, 62 = 0.22. 

Semi-quantitative agreement between experimental 
observation and the improved decay model allows us to 
conclude that this model reflects the basic physical mech¬ 
anisms responsible for the time dependence of the VLD. 


In particular, it accounts for the time delay in the deliv¬ 
ery of the energy flux from classical to quantum parts of 
superfluid turbulence. 


V. TOWARDS THEORY OF QUANTUM 
TURBULENCE WITH COUNTERFLOW 

In mechanically driven quantum turbulence, the mean 
velocities of the normal and superfluid components are 
known to coincide: Un = Ug. Numerous laboratory, nu¬ 
merical and analytical studies showed that under these 
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conditions the mutual friction between the normal and 
superfluid velocity components couple also their fluctua¬ 
tions: « Us{r,t) almost at all scales. This is not 

the case in thermally driven quantum turbulence, where 
the counterflow velocity Uns 7 ^ 0 partially decouples the 
normal and superfluid velocity fluctuations and enhances 
the turbulent energy dissipation due to the mutual fric¬ 
tion. In this Section we suggest a simple analytical model 
of the resulting energy balance in counterflow turbulence 
that predicts dramatic suppression of the energy spec¬ 
trum at intermediate and small scales. 

All ingredients of this model are already known. We 
just have to adopt Eqs. (13) for the one-dimensional en¬ 
ergy spectra Ss, Sn and to take in this equation expression 
for Ens{k) from Eq. (12). In this way we formulated the 
Differential Model for the normal and superfluid energy 
spectra Effk) and Es{k) of counterflow turbulence, given 
by the closed set of differential Eqs. (13), and (12). These 
equations may be solved numerically. 


which can be solved analytically with the boundary con¬ 
dition at (7 = 1 (i.e. k = kff)'. 



0 


Using now Eqs. (34b) and (35) one gets: 


(35a) 

(35b) 

(35c) 






(36b) 


A. Energy spectra in “symmetric” counterfiow 
turbnlence with pn = ps 


In order to get qualitative information on the form 
of the energy spectra Effk) and Es{k), we consider the 
derived equations in the important limiting case Elns 
7 ns, Then the steady-state version of Eqs. (13) 

and ( 12 ) simplify and take the form: 




(32) 


= a 




ViSs{k) + 


D{k)-Es,n} ■ 


Dramatic further simplification occurs for T ~ 1.95 K, 
when pn — Ps- In this case we expect Es(k) = Effk) = E 
and Eqs. (32) become an ordinary differential equation 
for E{k)-. 


=EtE{k) [D{k) - 1] , (33a) 

o CLfx 

Cq = q/Qx , Qx =2 El/{koVns) ■ (33b) 


Integral (35c) can be found analytically: 

4 


VStt 3 r 


— 2'\/3arctan['\/3 — 

-273arcta„[73 + 2,-7 + '“"'T"M 

zb/6 

3z^/^ j I 


- In 


1 - 


(I-f z2/3)2 
It has the following asymptotics: 


Simple analytical expression 
/(z)cs 


z4/3 

4 -b 0.7z2 


(36c) 


I{z) = 

2 ;^/^ / 62 ;^ \ 

for 2 ; <C 1, 

(37a) 

I{z) = 

VStt 3 

5 2z2/3 ’ 

for 2 : > 1 . 

(37b) 


(37c) 


approximates the exact Eq. (36c), with relative accuracy 
within 3% in the z < 2 region, while expression 


Here the dimensionless wave number q = k/ko with 
ko ~ tt/A and the outer scale of turbulence A is de¬ 
fined such that Eq. (32) are valid for q > 1. Introducing 
new dimensionless function 


'fiq) 

E{k) 


- 8^(fc/fco) i2 D? 


(34a) 

(34b) 


I{z) 


3 ^ TT 

“5 2 z 2/3 + 4 z 5/3 ’ 


(37d) 


works well for 2 > 2 . 

Notice that for ko < k kx the energy spectrum (36) 
deviates down from the K41 |-spectrum: 


we get from Eqs. (33a) the equation 
dq 


(34c) 


The spectrum for ki$> kx crucially depends on the value 
of A. There exists a critical value 

Acr = l//(oo) = 5 /( 7 rv^) « 0.92 . (38b) 
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For A < Acr the system asymptotically tends to K41 
spectrum 

<-) 

but with the energy flux = eo(l — A/Acr)^, smaller 
than the energy input rate £ 9 . The difference (eq — £ 00 ) 
is dissipated by mutual friction. This is similar to the 
subcritical LNV spectrum (18) of ^He turbulence with 
resting normal fluid component. 

For A > Acr^ £ik) = 0 for large k. In the differen¬ 
tial approximation used here, the spectrum E{k) sharply 
terminates at some finite in the same manner as the 
supercritical LNV ^He spectrum(19): 

The cutoff wave number /c* may be found from the equa¬ 
tion 

l = A[/(aj-/(Co)]«4l/(aj . (38e) 

When A —)• A^, —)■ 00 and £{k) cx k~^ at large k, 

exactly like in the critical LNV ^He spectrum (17). 

The energy spectra (36) for different values of A in 
the range 0 ^ ^ ^ 6 , that includes the critical value 
Acr ~ 0.92, are shown in Fig. 12a. We see that increase in 
the mutual friction force, characterized in Eq. (36) by the 
dimensionless parameter A, suppresses the energy spec¬ 
tra from the Kolmogorov-41 behavior £{k) oc k~^/^ (for 
A = 0 ) towards the critical spectrum £{k) oc k~^ at 
A = ^cr — 1- Further increase in A localizes energy 
spectra in the fc-space, as shown in Fig. 3(a). There is 
clear qualitative similarity of the energy spectra in con- 
terflow turbulence with the LNV spectra. This allowed 
us to use the LNV spectra in Sec. IV C for the analysis of 
the delay function Fdei(t) in our improved model of the 
VLB decay. 

B. Velocity structure functions in “symmetric” 
counterflow turbulence 

Recent visualization experiment in counterflow^^ re¬ 
ported that transversal velocity structure function 
'S' 2 ._l(?') oc r in an interval about one decade, adjacent 
to the outer scale of turbulence. To see how these obser¬ 
vations may be rationalized using our model, we consider 
the second order velocity structure functions, defined as 
follows: 

S 2 {r) = (|M(r) - m(O)P) , (39a) 

<S'2,_L(r') = ([^^(r) - u_l( 0)]^) , (39b) 

'S' 2 ,||(r') = ([u||(r) -u||(0)|2) . (39c) 

Here u± and M|| are projections of the turbulent velocity 
u on directions orthogonal and parallel to the separation 


r. In isotropic turbulence, assumed above in the analyt¬ 
ical model of the spectra, all structure functions depend 
only on r = |r| and all are proportional to each other. 
Up to a numerical factor of the order unity, they may be 
expressed via one-dimensional energy spectrum £(k) as 
follows: 

nOO 

52 (r) ~ / S{k)[l — s\n{kr)/(kr)]dk . (40a) 

Jo 

With the Kolmogorov-41 spectrum 5 k4i(A) = 

£o{ko/k)^^^ [Eq. (36) with A = 0], it gives the classical 
result 

52(r) cs 1.2U,2(^^^)2/3 ^ (40b) 

shown by the upper green dashed straight line in Fig. 12c. 
Next, we account for the fact that in reality the available 
range of k is limited: /cmin < k < fcmax- The value of 
fcmin ~ tt/A is determined by the outer scale of turbu¬ 
lence. For quantum turbulence in the superfluid com¬ 
ponent fcinax — In the Florida experiments^"^, our 

estimates show that fcmin — 0.3fcx and fcmax — 200fcx. 
Replacing limits in the integral (40a) by these values, we 
compute again S 2 (r) with the K41 spectrum- see upper 
green solid line in Fig. 12c. We observe the same scal¬ 
ing behavior S 2 {r) oc in the interval of about two 
decades (from r ~ 0.05 to r ~ 5). For nonzero values 
of A the log-log plots of S 2 (r) versus r can be consid¬ 
ered as approximately straight lines with the slope that 
increases with A. In particular, for A = 0.6 S 2 (r), shown 
in Fig. 12c by solid red line, it is practically indistinguish¬ 
able from the straight line with the slope -1-1 (shown by 
dashed red line) in the interval 0.02 < r < 2. This means 
that, for A = 0.6, S 2 {r) oc r with high accuracy in the 
interval of two decades. To see this better, in Fig. 12d we 
present the plot of S 2 {r) compensated by 1/r. The solid 
red line in Fig. 12(d) lays indeed very close to the black 
thin horizontal line. 

Notice that the energy spectrum for A = 0.6, used 
to find 52 (r) (red solid line in Fig. 12b), is essentially 
different from the scale-invariant spectrum £{k) oc k~‘^, 
shown in Fig. 12b by the dashed blue line. Using spec¬ 
trum £{k) oc we computed 52(r) again, with the 

result shown in Fig. 12d by the dashed blue line. Un¬ 
expectedly, this result demonstrates scale-invariant be¬ 
havior 52 (r) oc r on a shorter range. We also com¬ 
puted 52(r) using supercritical LNV spectrum (18), with 
fccr = l.lS/co- This spectrum, shown in Fig. 12b by the 
green dash-dotted line, is very different from the 1/fc^ 
behavior (blue dashed line). Nevertheless, the result¬ 
ing structure function 52 (r) (green dash-dotted line in 
Fig. 12d) again demonstrates the scale-invariant behav¬ 
ior of 52(r) oc r over more than two decades. 

We conclude that very different energy spectra, in¬ 
cluding the spectrum (36) with A ~ 0.6, found here, 
can result in the reported Ref. S 2 (r) oc r with some¬ 
what smaller extent of the scaling behavior, of about one 
decade. This means that our analytical model does not 
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Compensated by the K41 scaling superfluid energy spectra k^^^S{k) 


(a) 


(b) 
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FIG. 12: Color online. Panel (a) Log-Log plots of compensated normalized energy spectra A:®^^£i(A:)/£max for fco = 0 and with 
different values of A, starting with ^ = 0 (upper green line), through A — 0.2, 0.4, 0.6 (red thick line) 0.8,0.92 (thick cyan 
line - critical value), 1, 1.5, 2.0 3.0 and last line - j4 = 6 . Panel (b) Log-Log plots of compensated normalized energy spectrum 
fc®^®£’(fc)/£’niax for ko = O.Sfcx with A = 0.6 (solid red line), 1/k^ (after compensation by - blue dashed line and LNV 
subcritical spectrum fc“®[l -|- (fc/fc*)^^®]^ with fc* = l.lSfcx - green dash-dotted line. Panel (c) Log-log plots of the second order 
structure functions S 2 {r) (normalized by their large r limit), computed with Eq. (40a) and with the same values of A [and the 
same color code] as in Panel (a). The thick dashed straight lines indicate scaling laws: upper green line oc ), middle red 
line oc r and lower blue line oc r^). Panel(d) The normalized compensated by 1/r structure functions 52 (r), calculated using 
the spectra shown in Panel(b) (with the same color code). 


contradict the observation^’^. Nevertheless, at this stage 
we are not in the position to claim that the model ex¬ 
plains the observed simple behavior S 2 {r) (x. r. We will 
return to this point in conclusive Sec. VI. 

Here we notice only that, as the parameter A increases, 
the energy spectra become more and more localized at 
small k, as seen in Fig. 12(a), while the apparent slope of 
the corresponding structure functions increases, tending 
to 2. This is clearly seen in Fig. 12(c). The 52(r) for the 
largest H = 6 is shown by solid blue line in comparison 
with the scaling function r^, shown by the blue dashed 
line. The reason for such a behavior is simple: for large A, 
the energy spectrum terminates at small k and for scales 
smaller than 1/A:, the velocity field becomes smooth and 


differentiable. Thus the velocity field can be expanded in 
the Taylor series, the velocity difference across the sepa¬ 
ration r is proportional to r, and the structure function 
is proportional to r^. 

C. Approximations of the simple analytical model 

It should be stressed that in our approach we adopted 
some uncontrolled approximations and simplifications, 
widely used in the studies of classical hydrodynamic tur¬ 
bulence. Among them are celebrated hypotheses, sug¬ 
gested by Kolmogorov in 1941, concerning the small scale 
turbulent statistics: 
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1. universality (independence of the energy pumping); 

2. isotropy; 

3. locality of energy transfer over scale. 

These hypotheses have been justihed in numerous ex¬ 
periments, numerical and analytical studies of developed 
turbulence of classical fluids, mechanically driven at large 
scales. However, very little is known about turbulent 
statistics in the case of thermally driven quantum coun¬ 
terflow turbulence. 

We can quite easily accept the first hypothesis of the 
universality, for instance in the wide and long enough 
channels with reasonably well-controlled surface of the 
wall. 

The second assumption is the isotropy hypothesis. It is 
known from numerical simulations of (mainly space ho¬ 
mogeneous) counterflow turbulence that the anisotropy 
of the quantized vortex tangle is rather small (about 
10%), see e.g. Ref. However very little is 

known about anisotropy of the velocity fluctuations on 
scales R above the intervortex distance £. Without this 
knowledge and having clear understanding that the as¬ 
sumption of isotropy on scales R > i may be ques¬ 
tionable, we nevertheless assume isotropy of turbulence 
statistics from the very beginning. This simplifying 
assumption allows us to formulate a simple analytical 
model in terms of one-dimensional energy spectrum £{k), 
which is an angular average of the full (and possibly 
anisotropic) three-dimensional energy spectrum. 

The third assumption is the locality of the energy 
transfer, which is build-in in the algebraic closure (14b), 
used in our approach. Bearing in mind that the nonlinear 
terms in the two-fluid equations of motion for the normal 
and superfluid velocity components are the same as in 
the Navier-Stokes equation, and, moreover, if the result¬ 
ing energy spectrum is not very different from the clas¬ 
sical |-scaling for the classical fluid, we can accept also 
this assumption. The detailed analysis of the problem 
of locality in classical turbulence with the scale-invariant 
scaling £{k) oc shows'*®’^® that in the “window of lo¬ 
cality” 3 < a; < 5 the energy transfer is really local. For 
our case this means that for the subcritical spectra (with 
A < 1, when the the local slope x{k) = d\n[£{k)]/d\nk is 
within the window of locality), this assumption is reason¬ 
able. We also hope that the strong version of the locality 
assumption, the algebraic closure (14b), leads to quali¬ 
tatively correct energy spectra. The situation with the 
supercritical spectra is less simple. Definitely, the conse¬ 
quence of Eqs. (36) and (38d) that £{k) = 0 for k > k^, 
is an artefact of the algebraic closure. We think that for 
large k, the supercritical energy spectra will instead de¬ 
cay very fast with k, presumably oc k~^ with y > 5, as a 
result of the direct (non-local) energy transfer to the k- 
region from the intermediate region of scales about k-mt, 
where the local slope x{kint) — 5. 

Last but not least, an additional restriction of our ap¬ 
proach is full ignorance of possible long-living coherent 
structures at scales > £ that may contribute to the statis¬ 
tics of quantum turbulence with counterflow. There is a 


great deal of speculation in the studies of classical space- 
homogeneous hydrodynamic turbulence, but, neverthe¬ 
less, the question of their statistical relevance is still open. 
Moreover, the well justified multi-fractal models of clas¬ 
sical turbulence (see, e.g., textbook^^) describe in many 
details the statistics of classical turbulence without any 
direct reference to coherent structures. Coherent struc¬ 
tures in quantum turbulence can be characterized simply 
as “terra incognita”. 

In such a situation, much more experimental, numeri¬ 
cal and analytical work is required to formulate a theory 
of quantum turbulence which will account for the inter¬ 
play of coexisting classical and quantum forms of super¬ 
fluid turbulent energy in all relevant details. Neverthe¬ 
less, we consider our experimental findings and simple 
analytical models of steady-state and decaying quantum 
turbulence as a natural and perhaps even required step 
in a long way toward desired level of understanding and 
description of the basic physical mechanisms that govern 
quantum turbulence. 


VI. SUMMARY 

Being motivated by the challenge to understand the 
quantum turbulence occurring in superfluid ^He, we re¬ 
port this complementary experimental, numerical and 
theoretical treatise of turbulent coflow, counterflow and 
pure superflow of superfluid "^He in a channel. The level 
of agreement between the experimental observations and 
the analytical predictions for the time evolution of the 
vortex-line density in decaying turbulence, demonstrated 
in this paper, allows us to conclude that the developed 
basic and improved models adequately reflect the un¬ 
derlying physical processes responsible for the decay of 
quantum turbulence, originating from various types of 
steady superfluid "^He flows (coflow, counterflow, pure su¬ 
perflow), including 

- the interplay of classical and quantum processes, resul- 
sting in two decay laws of VLD; 

- the partial decoupling of the normal and superfluid 
velocity helds in pure superflow and counterflow turbu¬ 
lence, predicted in Ref. ; 

- the resulting suppression of energy spectra in these 
flows leading to the time-delay in the energy flux from 
classical to quantum length scales of turbulence in su- 
perffuid ^He. 

Being inspired by these findings, we made the hrst step 
towards the theory of steady-state, space homogeneous 
turbulence of counterflowing superfluid "^He. The sug¬ 
gested by us simple analytical theory results in the energy 
spectra £{k,A) given by Eq. (36) and shown in Fig. I2a. 
These spectra depend on the dimensionless parameter 
A, that describes the intensity of the mutual friction, re¬ 
sponsible for the coupling of the normal and superfluid 
velocity, relative to the counterflow velocity, responsible 
for their decoupling. For the particular value A ~ 0.6, 
the energy spectrum £{k,Q.6), shown in Fig. I2b is close 
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to 1/fc^ scaling and results in the second-order structure 
function 52 (r) oc r over the interval of about two decades, 
in agreement with the observations reported in Ref. 


APPENDIX A: SABRA-SHELL MODEL OF 
QUANTUM TURBULENCE 

A detailed study of the decay of large-scale turbulence 
has been done in the framework of the so-called Sabra- 
shell model of quantum turbulence 
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These equations represent a simplified version of the 
coarse-grained, two-fluid, gradually-damped Hall-Vinen- 
Bekarevich-Khalatnikov (HVBK) equations in the k- 
representation. They mimic the statistical behavior of 
fc-Fourier components of the turbulent superfluid and 
normal velocity fields in the entire shell of wave vectors 
km < k < km+i by complex shell velocity The 

shell wave numbers are chosen as a geometric progres¬ 
sion km = fcoA™, where m = 1,2,... M are the shell 
indexes, and we have used the shell-spacing parameter 
A = 2, /co = 1/16 and M = 28 shells. 

Similarly to the HVBK (and the Navier-Stokes or Euler 
equation), the NL[um] term in Eq. (Ale) is quadratic in 
velocities, proportional to k and conserves (in the force¬ 
less, inviscid limit) the kinetic energy E = \ J^m 
provided that a -|- & -I- c = 0. We used here the Sabra ver¬ 
sion^® of NL[Mm] with the traditional (and physically mo¬ 
tivated) choice b = c = —a/2, which describes important 
features of superfluid turbulence including intermittency 
corrections 

In Eq. (Ala), Vn = ^/Pn is the kinematic viscosity of 
the normal component, i.e., its dynamical viscosity /i, 
normalized by the normal fluid density. The effective 
superfluid viscosity t's describes the energy sink in super¬ 
fluids, e.g. due to the vortex reconnections. For more 
details of the origin and role of see Ref. The mu¬ 
tual friction terms F)))’®, given by Eq (Aid), is just the 
fc-Fourier transform of Eqs. (10a) and (10b). 

As a first step in our study of the large-scale turbulence 
decay we can simplify the procedure further, using the 
fact that the turnover time in our situations is longer than 
the coupling time 1/Qns given by Eq. (lOd). Therefore 
the first stage of the decay, during which the normal- 
and superfluid velocities become coupled, is short with 
respect to the time required for the developing of the 


Richardson-Kolmogorov cascade. Skipping this stage, we 
can assume that the velocities are fully coupled according 
to Eq. (8). In this case Eqs. (Al) turn into one Sabra- 
equation for Um = Um = 


d ,2 


= NL[um], 


(A2) 


with V = [jiA- i'sPs)/p- 

The equations (A2) were solved using the adaptive 
time step 4th order Runge-Kutta with exponential time 
differencing^®. For more details see Ref. The evo¬ 
lution of the system was followed for about 10 ®ta. In all 
cases we perform 10^ simulations with the same initial 
total energy, but different, randomly distributed phases 
of initial shell-velocities and perform ensemble averaging 
over initial conditions. 


APPENDIX B: SOME DEFINITIONS AND 
KNOWN RELATIONSHIPS 

To define the one-dimensional energy spectra En(k), 
Es{k) and cross-correlation 5ns(fc) we need to recall some 
definitions and relationships, that are well-known in sta¬ 
tistical physics. The first is the set of Fourier transforms 
in the following normalization: 


K,sif',t) = 

[ (27r)3 '*^n,s(fc,t)exp(ife • r), 

(Bla) 

III 
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/ ^ exp{-iujt), 

(Bib) 

Vn,s{k,Uj) = 

j drdt Mn_s(r, t) exp[i(a;t — fe • r)] 

. (Blc) 


The same normalization is used for other objects of in¬ 
terest. 

Next we define the simultaneous correlations and cross¬ 
correlations in fe-representation, (proportional to 5{k — 
k') due to homogeneity): 

(Un(fc,t)-<(fc',t)) = (2^)®Fn„(fe)d(fc-fc'),(B2a) 

(u,(fc,t)-<(fc',t)) = (2^)®F,,(fc)(5(fc-fc'), (B2b) 
(un(fc,t)-<(fc',t)) = (2^)®Fn,(fe)d(fe-fc') • (B2c) 

We also need to define cross-correlations {v^ ■ v*) in 
(fe, u;)-representation: 

{Vn{k,uj) ■ V*{k',uj')) (B3a) 

= (27r)‘‘F„s(fe,a;)d(fe-fe')<5(w-a;') . 

This object is related to the simultaneous {v^ ■ v* ) cross¬ 
correlation (B2c) via the frequency integral: 

{Vn{k,t) ■ v*{k',t)) = J dwFns(fe,w) • (B3b) 

Here and below “tilde” marks the objects defined in 
(fe, a;)-representation. 
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It is known also that the fc-integration of the correla¬ 
tions (B2) produces their one-point second moment: 
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^Ess{k,t) 

= (l«s(r',t)|3) , 

(B4b) 

-Ens{k,t) 

= {un{r,t) ■ Us{r,t)) . 

, (B4c) 
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